Final spins from the merger of precessing binary black holes 
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The inspiral of binary black holes is governed by gravitational radiation reaction at binary sepa- 
rations r < 1000M, yet it is too computationally expensive to begin numerical-relativity simulations 
with initial separations r > 10M. Fortunately, binary evolution between these separations is well 
described by post-Newtonian equations of motion. We examine how this post-Newtonian evolu- 
tion affects the distribution of spin orientations at separations r ~ 10M where numerical-relativity 
simulations typically begin. Although isotropic spin distributions at r ~ lOOOAf remain isotropic 
at r ~ 10M, distributions that are initially partially aligned with the orbital angular momentum 
can be significantly distorted during the post-Newtonian inspiral. Spin precession tends to align 
(anti-align) the binary black hole spins with each other if the spin of the more massive black hole 
is initially partially aligned (anti-aligned) with the orbital angular momentum, thus increasing (de- 
creasing) the average final spin. Spin precession is stronger for comparable-mass binaries, and could 
produce significant spin alignment before merger for both supermassive and stellar-mass black hole 
binaries. We also point out that precession induces an intrinsic accuracy limitation (< 0.03 in the 
dimensionless spin magnitude, < 20° in the direction) in predicting the final spin resulting from the 
merger of widely separated binaries. 

PACS numbers: 04.25. dg, 04.25.Nx, 04.70.-s, 04.30.Tv 



I. INTRODUCTION 

The existence of black holes is a fundamental pre- 
diction of general relativity. Isolated individual black 
holes are stationary solutions to Einstein's equations, 
but binary black holes (BBHs) can inspiral and even- 
tually merge. BBH mergers offer a unique opportunity 
to test general relativity in the strong- field limit, and 
as such are a primary science target for current and fu- 
ture gravitational-wave (GW) observatories like LIGO, 
VIRGO, LISA, and the Einstein telescope. BBH merg- 
ers are also important for cosmology, as they can serve 
as standard candles to help determine the geometry and 
hence energy content of the universe [l], Q ■ Astrophysi- 
cal BBHs are found on at least two very different mass 
scales. Compact objects believed to be stellar-mass black 
holes have been observed in binary systems with more 
luminous companions. These black holes are the rem- 
nants of massive main-sequence stars, and binary systems 
with two such stars may ultimately evolve into BBHs. 
On larger scales, supermassive black holes (SBHs) with 
masses 10 6 < M/M Q < 10 9 reside in the centers of most 
galaxies. They can be observed through their dynamical 
influence on surrounding gas and stars, and when accret- 
ing as active galactic nuclei (AGN). SBHs will form bi- 
naries as well, following the merger of two galaxies which 



each host an SBH at their center. 

In order to merge, BBHs must find a way to shed their 
orbital angular momentum. At large separations, bi- 
nary SBHs will be escorted inwards by dynamical fric- 
tion between their host galaxies The BBHs be- 
come gravitationally bound when the sum of their masses 
AI = mi +rri2 exceeds the mass of gas and stars enclosed 
by their orbit. The binary hardens further by scattering 
stars on "loss-cone" orbits that pass within a critical ra- 
dius [If, though this scattering may stall at separations 
r ~ 0.01 — 1 pc unless these orbits are refilled by stellar 
diffusion [f|. Unlike stars, gas can cool to form a cir- 
cumbinary disk about the BBHs. A circumbinary disk of 
mass Md and radius will exert a tidal torque 
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on the binary in the limit that the BBH mass ratio 
q = milm\ < 1 is small and | — r| <r Through- 
out this paper we use relativists' units in which New- 
ton's constant G and the speed of light c are unity. At a 
sufficiently small separation rewj the magnitude of this 
tidal torque will fall below that of the radiation-reaction 
torque [9| 



w 



32t? 2 M 9 / 2 
5r 7 /2 



(1.2) 



"Electronic address: kcsdcn@tapir.caltech.edu 



^Electronic address: spcrhake@tapir.caltcch.edu 
t Electronic address: berti@phy.olemiss.edu 



where r\ = mxra^jM' 1 is the symmetric mass ratio. Once 
Tqw > Td, the inspiral of the BBH is dominated by 
radiation reaction. The precise value of row depends on 
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the properties of the circumbinary disk, but an order-of- 
magnitude estimate is given by 
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validity of the PN expansion for spinning precessing bi- 
naries, see e.g. [HI HH)- Numerical relativists can sim- 



where Mg is the mass of the larger black hole in units of 
10 8 M©, th is the dynamical friction timescale for a hard 
binary, and t gas is the evolution timescale from gaseous 
tidal torques. 

General relativity completely determines the inspiral 
of BBH systems from separations less than tqw- These 
systems are fully specified by 7 parameters: the mass ra- 
tio q and the 3 components of each dimensionless spin 
Xi 2 = Si, 2/1^1 2- To a good approximation the indi- 
vidual masses and spin magnitudes xi,2 = \Xi 2 1 remain 
constant during the inspiral, so only the precession of 
the two spin directions needs to be calculated. At an 
initial separation r% = 1000 M ~ row, the binary's or- 
bital speed v/c <C 1 and the spin-precession equations 
can therefore be expanded in this small post-Newtonian 
(PN) parameter. The PN expansion remains valid until 
the BBHs reach a final separation rj = 10M, after which 
their evolution can only be described by fully nonlinear 
numerical relativity (for more precise assessments of the 
3 N expa 
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ulate BBH mergers from separations tnr — rj [12Hl^ 
but these simulations are too computationally expensive 
to begin when the binaries are much more widely sepa- 
rated. The GWs produced in the merger and the mass, 
spin, and recoil velocity of the final black hole depend 
sensitively on the orientation of the BBH spins at Tnr, 
so it is important to determine what BBH spin orienta- 
tions are expected at ri and whether these orientations 
are modified by the PN evolution between and r/. 

The answer to the first of these questions comes from 
astrophysics, not general relativity. At very large sepa- 
rations, the two black holes are unaffected by each other 
and one would therefore expect an isotropic distribution 
of spin directions. However, an isotropic distribution of 
spins at rj would imply that most mergers would re- 
sult in a gravitational recoil of ~ 1000 km/s for the final 
black hole Recoils this large would eject SBHs 

from all but the most massive host galaxies [l8[ , in seem- 
ing contradiction to the observed tight correlations be- 
tween SBHs and their hosts [l9l-l2l| . This problem can be 
avoided if Lcnsc-Thirring precession and viscous torques 
align the spins of the BBHs with the accretion disk re- 
sponsible for their inwards migration [2214241 ]. The effi- 
ciency of this alignment depends on the properties of the 
accretion disk, but -/V-body simulations using smoothed- 
particle hydrodynamics (SPH) suggest that the residual 
misalignment of the BBH spins with their accretion disk 
at r, could typically be ~ 10° (30°) for cold (hot) accre- 
tion disks 17 1 . 



The second question, does the distribution of spin di- 
rections change as the BBHs inspiral from j-j to r/, can 
be answered by evolving this distribution over this in- 
terval using the PN spin-precession equations. We will 
describe these PN equations and our numerical solutions 
to them in Sec. [H] The precession of a given spin configu- 
ration in the PN regime can be understood in terms of the 
proximity of that confi gura tion to the nearest spin-orbit 
resonance. Schnittman [25j identified a set of equilibrium 
spin configurations in which both black hole spins and the 
orbital angular momentum lie in a plane, along with the 
total angular momentum J = L + m\xi + m 2%2- Ti the 
absence of radiation reaction, J is conserved. For these 
equilibrium configurations, the spins and orbital angular 
momentum remain coplanar and precess jointly about J 
with the angles $1,2 between L and Xi 2 remaining fixed. 
The equilibrium configurations can thus be understood 
as spin-orbit resonances since the precession frequencies 
of L and Xi 2 about J are all the same. Once radiation re- 
action is added, the spins and orbital angular momentum 
remain coplanar as the BBHs inspiral, although Q\ and 
6*2 slowly change on the inspiral timescale. Not only do 
resonant configurations remain resonant, but configura- 
tions near resonance can be captured into resonance dur- 
ing the inspiral. The resonances are thus very important 
for understanding the evolution of generic BBH systems, 
although the resonances themselves only occupy a small 
portion of the 7-dimensional parameter space character- 
izing generic mergers. We shall review these spin-orbit 
resonances in more detail in Sec. IIIII 

Bogdanovic et al. [23| briefly considered whether spin- 
orbit resonances could effectively align SBH spins with 
the orbital angular momentum following the merger of 
gas-poor galaxies. They found that for a mass ratio 
q = 9/11 and maximal spins Xi = X2 — 1, an isotropic 
distribution of spins at 7"j = 1000M remains isotropically 
distributed when evolved to 77 = 10M. They there- 
fore concluded that an alternative mechanism, such as 
the accretion torques considered later in their paper, is 
needed to align the BBH spins with L. This conclu- 
sion is supported by a much larger set of PN inspirals 
presented by Herrmann et al. [26j who found that for 
equal- mass BBHs, an isotropic distribution of spins at 
40M yields a flat distribution in 0036*12 at 7AM. Here 
and in this paper #12 is the angle between the two spins 
Xi and v 9 . In the final plot of their paper, Herrmann 
et al. [261 ] revealed their discovery of an anti-correlation 
between the initial and final values of cos 6*12 for q = 2/3 
BBHs with equal dimensionless spins xi — Xi = 0.05. 
Investigation of this anti-correlation was left to future 
work. Lousto et al. [27| also found indications that 
an initially isotropic distribution of spins can become 
non-isotropic during the PN stage of the inspiral. For 
a range of mass ratios l/16<g<l and equal spins 
Xi = X2 = (0.485,0.686,0.97), they found that an 
isotropic spin distribution at 50M develops a slight but 
statistically significant tendency towards anti-alignment 
with the orbital angular momentum L. This amplitude 
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of anti-alignment scales linearly in the BBH spin magni- 
tudes and appears to decrease as q — > 0. 

We perform our own study of PN spin evolution from 
Ti to Tf for several reasons. BBHs get locked into spin- 
orbit resonances at a separation 
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which can become large in the equal- mass (q — > 1) limit 
[25| . This limit is important, as the largest recoil ve- 
locities occur for nearly equal-mass mergers. Numerical 
integration of the PN equations has shown that for a 
mass ratio q = 9/11, spin-orbit resonances affect spin 
orientations at separations r ~ 100071/. This is a much 
larg er separation than was considered in previous studies 
[26l l27| of spin alignment, which may therefore have failed 
to capture the full magnitude of the effect. These studies 
also focused on whether an initially isotropic distribution 
of spins becomes anisotropic just prior to merger. How- 
ever, as discussed above, tidal torques from a circumbi- 
nary disk partially align spins with the orbital angular 
momentum at separations r ?*gw before relativistic ef- 
fects become important. As we will show in Sec. lIVi such 
partially aligned distributions can be strongly affected by 
spin-orbit resonances despite the fact that isotropic dis- 
tributions remain nearly isotropic. We will consider how 
spin precession affects the final spin magnitudes and di- 
rections in Sec. |Vj The evolution of the distribution of 
BBH spin directions between and r/ changes the dis- 
tribution of final spin magnitudes and directions from 
what it would have been in the absence of precession. In 
addition, spin precession introduces a fundamental un- 
certainty in predicting the final spin of a given BBH sys- 
tem. At large separations, a small uncertainty in the 
separation leads to an uncertainty in the predicted time 
until merger that exceeds the precession time. In this 
case, one cannot predict at what phase of the spin pre- 
cession the merger will occur and thus the resulting final 
spin. We will explore this uncertainty in Sec. lVIl A brief 
discussion of the chief findings of this paper is given in 
Sec. HTTP 



der, and spin effects are included up to 2PN order. These 
evolution equations were chosen for consistency with pre- 
vious workjin particular with the study by Barausse and 
Rezzolla [30( of the final spin resulting from the coa- 
lescence of BBHs and with the statistical investigation 
of spinning BBH evolutions using Graphics Processing 
Units by Herrmann et al. prjj . Lousto et al. 27 1 evolved 
a large sample of spinning BBH systems using a non- 
rcsummed, PN expanded Hamiltonian. The convergence 
properties of non-resummed Hamiltonians for spinning 
BBH systems are somewhat problematic (see e.g. Fig. 1 
of Ref. [l(| ) , and it will be interesting to repeat these sta- 
tistical investigations of precessing BBH systems using 
the effective-one-body resummations of the PN Hamilto- 
nian recently proposed by Barausse et al. jU, HH • 



In our simulations, the spins evolve according to 
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are the spin precession frequencies averaged over a circu- 
lar orbit, including the quadrupolc-monopolc interaction 
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is the Newtonian orbital angular momentum, and 
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II. POST-NEWTONIAN EVOLUTION 

We evolve spinning BBH systems along a sequence of 
quasi-circular orbits according to the PN equations of 
motion for precessing binaries first derived by Kidder 
HU, and later used by Buonanno, Chen and Vallisncri to 
build matched-filtering template families for GW detec- 
tion [2^] . The adiabatic evolution of the binary's orbital 
frequency is described including terms up to 3.5PN or- 



is the orbital frequency. In the absence of gravitational 
radiation, J and |Ljv| are constant, implying that the 
direction of the orbital angular momentum evolves ac- 
cording to 
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where S = Si + S 2 . Once radiation reaction is included, 
the orbital frequency slowly evolves as 
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where je — 0.577 is Eulcr's constant, 6 
and we have defined 
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The two terms in square parentheses on the third line 
of Eq. (1 2.61) are due to the quadrupole-monopole interac- 
tion j34| and to the spin-spin self interaction [35| , respec- 
tively, and they were neglected in the statistical study of 
Ref. [26|. Their sum agrees with Eq. (5.17) of Ref. [361 ) . 

The numerical integration of this system of ordinary 
differential equations is performed using the adaptive 
stepsize integrator StepperDopr5 J37|. The evolution 
of any given BBH system is specified by the following 
parameters: the initial orbital frequency u>i, the binary's 
mass ratio q = mz/rn-L, the dimensionless magnitude of 
each spin \ii and the relative orientation (9i, </>j) of each 
spin with respect to the orbital angular momentum at 
time t = (i = 1, 2). To monitor the variables along 
the whole evolution we output all quantities using a con- 
stant logarithmic spacing in the orbital frequency at low 
frequencies, and the stepsize as used in the integrator 
at high frequencies. Typically this results in a total of 
about 64,000 points in the range Muj € [Mu>i, Mcjf], 
where Mu, = 3.16 x 10" 5 and Mw/ = 0.1. Numeri- 
cal experimentation indicates that a tolerance parameter 
ATOL= 2 x 10~ 8 in the adaptive stepsize integrator is suf- 
ficient for a pointwise accuracy of order 1% or better in 
the final quantities. Therefore the error induced by the 
numerical integrations of the PN equations of motion is 
subdominant with respect to the errors induced by pre- 
cessional effects and by fits of the numerical simulations, 
which will be one of the main topics of this paper. 



III. SPIN-ORBIT RESONANCES 

In this Section, we review the equilibrium configura- 
tions of BBH spins first presented in Schnittman [25j for 



which the Newtonian orbital angular momentum L^r and 
individual spins Si, 2 all precess at the same resonant 
frequency. As discussed briefly in the Introduction, at 
a given binary separation r fully general quasi-circular 
BBHs are described by 7 parameters: the mass ratio q 
and the 3 components of each black hole spin. In spheri- 
cal coordinates with Ljv defining the z-axis, each spin is 
given by its magnitude Si = m 2 Xi and direction (9i, (j>i) 
(i = 1, 2). In the PN limit for which this analysis is valid, 
a clear hierarchy 
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exists between the orbital time i or b oc r 3 / 2 , the pre- 
cession time t p ~ r 5 / 2 , and the radiation time 

^gw ~ Eqvj/E cx r 4 . This hierarchy implies that the 
BBH spins will precess many times before merger leav- 
ing only their relative angular separation A(f> = <p2 — 4>i 
in the orbital plane well defined. This reduces the BBH 
parameter space to 6 dimensions. Since the mass ratio 
and individual spin magnitudes are preserved during the 
inspiral, a given BBH evolves through the 3-dimensional 
parameter space {6\ , 2 , A</>) on the precession timescale 
t p . This evolution is governed by the spin precession 
equations (|2.ip . 

Schnittman (25| discovered a one-parameter family 
of equilibrium solutions to these equations for which 
(0i,02, Ac/)) remain fixed on the precession timescale t p . 
These solutions have A</> = 0° or 180°, implying that 
Ljv, Si and S2 all lie in a plane and precess at the same 
resonant frequency about the total angular momentum 
J, which remains fixed in the absence of gravitational ra- 
diation. The values of 6\ t 2 for these resonances can be 
determined by requiring the first and second time deriva- 
tives of Si • S2 to vanish. This is equivalent to satisfying 
the algebraic constraint 

(fii x Si) • [S 2 x (Ljv + SO] 

= (fi a x S 2 ) • [Si x (Ljv + S 2 )] . (3.2) 



Since Ljv appears in Eq. (|3.2p both explicitly and implic- 
itly through f2i. 2 , the resonant values of t9i )2 vary with 
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FIG. 1: Spin-orbit resonances for maximally spinning BBHs 
with a mass ratio of q — 9/11. The dotted black diag- 
onal indicates where 8i = 82- Solid black curves below 
(above) this diagonal show (61, 62) for the one-parameter fam- 
ilies of equilibrium spin configurations with A<f> = 0°(180°) 
at different fixed binary separations. Approaching the di- 
agonal from below, these curves correspond to separations 
r = 1000, 500, 250, 100, 50, 10M. The curves approaching 
from above correspond to separations r = 250, 50, 20, 10M. 
The long-dashed red curves show how #1,2 evolve as mem- 
bers of these resonant families inspiral from = 1000M to 
tj = 10M. The projection S ■ Lat of the total spin S onto the 
orbital angular momentum Lat is constant along the short- 
dashed blue lines, while the projection So • Ljv of the EOB 
spin So is constant along the dot-dashed green lines. 



the binary separation. This is crucial, as otherwise these 
one-parameter families of resonances would affect only 
a small portion of the 3-dimensional parameter space 
(61 , 02 , A<^>) through which generic BBH configurations 
evolve. As gravitational radiation slowly extracts angular 
momentum from the binary on the radiation time tew, 
the resonances sweep through a significant portion of the 
#2) plane. The angular separation A0 of a generic 
BBH is varying on the much shorter precession time t p , 
and thus has a significant chance to closely approach the 
resonant values A<j> = 0° or 180° at some point during 
the long inspiral. Such generic BBHs will be strongly in- 
fluenced or even captured by the spin-orbit resonances, 
as we will see in detail in Sec. HVl 



We show the dependence of the spin-orbit resonances 
on r for maximally spinning BBHs in Figs.Uand El T hose 
resonances with A</> = 0° (shown in Fig. 2 of [2q|) al- 
ways have 6\ < 02, and thus appear below the diagonal 



FIG. 2: Spin-orbit resonances for maximally spinning BBHs 
with a mass ratio of q = 1/3. Other than the different mass 
ratio, this figure is very similar to Fig. [T] The solid black 
curves approaching the diagonal from below correspond to the 
families of resonant spin configurations at r = 50, 20, 10, 5M, 
while those approaching from above correspond to separations 
r = 20, 10, 5M. 



cos 0i = cos 6*2 in our Figs. Q] and [21 Those resonances 
with A<j) = 180° (shown in Fig. 3 of HH) have 6>i > 2 
and therefore appear above the diagonal in our Figs. [T] 
and [2] We plot (cos 9±, cos 9 2 ) rather than {61,62) like 
PH because isotropically oriented spins should have a 
flat distribution in these variables. 

In the limit r — > 00, so that also |Ljv| — > 00, the reso- 
nant configurations either have Si or S2 aligned or anti- 
aligned with (either 6\ or 62 equals to 0° or 180°). 
This corresponds to the four edges of the plot in Fig. [T] 
For smaller fixed values of |Ljv|, the values {61,62) for 
the one-parameter families of resonant configurations ap- 
proach the diagonal 6\ = 62- BBHs in spin-orbit reso- 
nances at large values of |Ljv| (large r) remain resonant 
as they inspiral. As gravitational radiation carries away 
angular momentum, r decreases and 0± t 2 for individual 
resonant BBHs evolves towards this diagonal along the 
red long-dashed curves in Fig. [1] For resonances with 
A(/> = 0° (those below the diagonal), this evolution aligns 
the two spins with each other. Symmetry implies that 
aligning the spins with each other will lead to larger final 
spins and smaller recoil velocities [38|, H^] . 

The projection 

S • Ljv = Si cos 6>i + S 2 cos6 2 (3.3) 
of the total spin S = Si + S2 parallel to the orbital an- 
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gular momentum is constant along the short-dashed blue 
lines in Figs. Q] and O These blue lines have steeper slopes 
than the red lines along which the resonant binaries in- 
spiral. This implies that the total spin S becomes anti- 
aligned (aligned) with the orbital angular momentum for 
resonant configurations with Acf> = 0°(180°), leading to 
smaller (larger) final spins. The families of resonances 
with A<j> = 0° (below the diagonal) sweep through a 
larger area of the (cos 0\, cos 62) plane as the BBHs in- 
spiral, and approach the diagonal more closely. This 
implies that anti-alignment may be more effective than 
alignment, which might explain the "small but statisti- 
cally significant bias of the distribution towards counter- 
alignment" in S-Ljv noted in Lousto et al. [27| . However, 
Table IV of [13] indicates that both Si and S2 individu- 
ally become anti-aligned with L^v , whereas the spin-orbit 
resonances would align one black hole while anti-alig ning 



the other. All of the PN evolutions in Lousto et al. [27 1 
begin at separations of r = 50 M, which corresponds to 
the A(j> = 0° curve in Fig. [T] that is second closest to 
the diagonal. The resonances sweep through most of the 
plane below the diagonal at larger separations, suggesting 
that these short-duration PN evolutions may have failed 
to capture the full magnitude of the anti-alignment. We 
will investigate this possibility in Sec. IIVI 

Another interesting feature of Figs. Q] and [2] is that the 
red long-dashed curves along which the BBHs inspiral are 
nearly parallel to the dot-dashed green lines along which 
the projection So • Ljv of the effective-one-body (EOB) 
spin 



S = (l + <z)S 1 + (l + g- 1 )S 2 



(3.4) 



is constant. The conservation of this quantity at 2PN 
order was first noted in Ref. [33| and follows directly from 
Eqs. ([23]). (|2~2|) . and (|23)l . The conservation of S ■ L N 
rather than S ■ L jv itself allows for the possible alignment 
of the total spin S discussed in the previous paragraph. 

We conclude this Section by briefly discussing how the 
spin-orbit resonances vary with the mass ratio q, as can 
be seen by comparing the q = 9/11 resonances in Fig. [1] 
with the (7=1/3 resonances in Fig. [2j The most pro- 
nounced differences are that the q = 1/3 resonances 
sweep away from the edges of the (cos #1, cos #2) plane 
at much smaller values of the separation r, and do not 
approach the diagonal as closely. This is consistent with 
the decreasing value of ri oc k in Eq. (|1.4p as q — > 0. In this 
limit both t p and tew are proportional to q^ 1 , implying 
that generic BBHs will be less likely to be affected by 
the resonances as they sweep through the plane over a 
smaller range in r. BBHs already in a resonant configu- 
ration will also be less affected since the resonant curves 
do not approach the diagonal as closely. The red long- 
dashed curves showing the inspiral of resonant configu- 
rations have steeper slopes for q = 1/3, consistent with 
the larger black hole being immune to its smaller com- 
panion in the limit q — > 0. This seems to contradict the 
puzzling result presented in Table IV of Lousto et al. [27} 
that it is the smaller companion that remains randomly 
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FIG. 3: Distributions of (cos Si, cos 82) at different separations 
r for 1000 initially isotropic maximally spinning BBHs with 
a mass ratio q — 9/11. The top left panel shows the initial 
10 x 10 x 10 grid, evenly spaced in (cos 81, cos 82, At/)). The 
dotted vertical lines show cos 81 — ±0.4. The 300 blue squares 
initially have cos 8± > 0.4, the 400 green triangles initially 
have —0.4 < cos Si < 0.4, and the 300 red circles initially 
have cos Si < —0.4. The values of (81,82) for these BBHs are 
shown in the top right, bottom left, and bottom right panels 
after they have inspiraled to separations of r = 1000, 100 and 
10M respectively. 



distributed during the inspiral. We will examine this be- 
havior as well in the next Section. 



IV. SPIN ALIGNMENT 

In this Section, we examine the extent to which the 
spins of generic (i.e. misaligned) BBH configurations be- 
come aligned with the orbital angular momentum and 
each other as the BBHs inspiral from r, = 1000M to 
r/ = 1QAI. Although we use maximally spinning BBHs 
to demonstrate this alignment, the magnitude of the 
alignment is comparable for all BBHs with \i,2 ^ 0.5 
as shown in Fig. 11 of [25|. We first consider initial 
spin configurations given by a uniform 10 x 10 x 10 grid 
evenly spaced in (cos£q, cos #2, A<p). This distribution is 
isotropic, and would be expected in the absence of an 
astrophysical mechanism to align the spins. BBHs with 
isotropically oriented spins might form in gas-poor merg- 
ers of SBHs and mergers of stellar-mass black holes in 
dense clusters. 

In Fig. [31 we show how the distribution of 
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FIG. 4: Distributions of (cos #i, cos 62) at different separations 
r for 1000 initially isotropic maximally spinning BBHs with 
a mass ratio q — 1/3. The different panels, points, and lines 
are the same as those given for q — 9/11 in Fig. [3] 



(cos 9i, cos 62 ) evolves as maximally spinning BBHs with 
a mass ratio q = 9/11 inspiral from slightly beyond 
n = 1000 M tor/ = 10M. The top left panel shows 
our initial evenly spaced 10 x 10 x 10 grid. The points 
are colored to indicate their initial value of cos 81 : blue 
squares begin with cos#i > 0.4 (#1 < 66°), green tri- 
angles with —0.4 < cos 6*i < 0.4, and red circles with 
cos 0i < —0.4. The dotted vertical lines cos# = ±0.4 
denote these boundaries. Only 100 points are visible in 
the top left panel, as the different values of A<fi cannot 
be distinguished in this two-dimensional projection. Spin 
precession reveals all 1000 points after the BBHs have in- 
spiralcd to n = 1000M as seen in the top right panel. 
Notice that the spins of all 1000 BBHs precess in a way 
that conserves the projection of So onto (parallel 
to the dot-dashed green lines in Fig. [1]). This is not a 
special feature of the spin-orbit resonances, but occurs 
for generically oriented spins as well. These generic spin 
configurations do not individually preserve (cos 0\, cos 62 ) 
over a precession time t p like the resonant configurations 
do, but they do preserve the combination So ■ Ljv- This 
precession continues as the BBHs inspiral to r = 100A/ 
and Tf = 10M as shown in the bottom left and bot- 
tom right panels of Fig. [3] By the time they reach 
rf = 10M the green points have diffused to fill most 
of the (cos 0i, cos 9 2 ) plane, while the blue (red) points 
have diffused into the upper right (lower left) portion 
of the middle —0.4 < cos#i < 0.4 region. The bottom 



right panel, if the points had not been colored, would re- 
produce Fig. 1 of Bogdanovic et al. HH and therefore 
support their conclusion that isotropically distributed 
spins remain isotropic as they inspiral. However, the col- 
ors reveal that PN evolution can drastically alter spin 
distributions that have been partially aligned by a cir- 
cumbinary disk. For example, if the spin of the more 
massive black hole was aligned so that cos#i > 0.4 at 
ri = 1000M (shown by our blue points), by the time the 
binary reached r/ = 10M the larger spin could easily lie 
in the orbital plane and thus give rise to a smaller final 
spin and potentially large "superkick" [l5|, . 

For comparison, we show the inspiral of the same 
10 x 10 x 10 grid of maximally spinning BBHs with a 
mass ratio q = 1/3 in Fig. SJ The points diffuse along 
the steeper lines that preserve So • L jv for this less equal 
mass ratio. This inhibits their ability to diffuse across 
the cos B\ = ±0.4 boundaries, again shown by the ver- 
tical dotted lines. Even at r/ = 10M only a few points 
have trickled between the three regions. Since the spin of 
the more massive black hole remains aligned with the or- 
bital angular momentum, one would expect a large final 
spin and an absence of superkicks for such small mass 
ratios. We will examine in detail how spin alignment 
affects recoil-velocity distributions in future work. 

In Fig. [5] we show how the joint probability distribu- 
tion function for Acj) and cos 0\2 evolves for our evenly 
spaced 10 x 10 x 10 grid of initially isotropic BBH spin 
configurations. As defined in the Introduction, cos #12 is 
the cosine of the angle between Si and S2. It can be 
expressed in terms of the individual spin angles as 



cos 6*12 = sin 6\ sin 82 cos A</> + cos Q\ cos 6*2 



(4.1) 



and has a flat distribution between -1 and 1 for isotropic, 
uncorrelated spins such as those given by our 10x10x10 
grid. However, as seen in Eq. (|4.1[) . the values of cos Q\2 
and cosA</> are correlated; for a given value of Acj> the 
distribution of cos 6*12 is peaked about cosA</> for flat 
distributions of cos 9\ and cos 82 . This can be seen in 
Fig. [5] from the clustering of points about the curve 
cos 6*12 — cos Acj). Although cos 6*12 and cos A(j> are corre- 
lated even for isotropic spins, geometry implies that both 
are initially uncorrelated with the value of cos 6\ . This is 
revealed by the identical distributions of the red, green, 
and blue points in the top left panel of Fig.[5]to within the 
resolution of our grid. These distributions do not remain 
identical as the BBHs inspiral from r.; = 1000M to rj = 
10M. Influenced by the A<fi = 0° spin-orbit resonances 
below the diagonal in Fig. [TJ the blue points become con- 
centrated about A<f> = 0°, cos 6*12 = 1 by the time they 
reach rj. The red points, similarly influenced by the 
Atfi = ±180° resonances above the diagonal in Fig. [TJ 
become concentrated about A<j) = ±180°, cos 6>i 2 = — 1. 
The effect of this spin alignment on the spin of the final 
black hole will be explored in detail in the next Section, 
while the effect on recoil velocities will be examined in fu- 
ture work. Qualitatively, alignment of the spins with each 
other (cos 0i2 —5- 1) increases the final spin and reduces 
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FIG. 5: Distributions of (A<f>, cos 612) at different separations 
r for 1000 initially isotropic maximally spinning BBHs with 
a mass ratio q — 9/11. The top left panel shows the initial 
10 x 10 x 10 grid of BBH spin configurations, evenly spaced 
in (cos 61, cos 82, This distribution is peaked about the 

curve cos 612 = cos A<f) shown by the dotted curve. The points 
are colored according to their initial values of cos 0i as in 
Fig. [3] The top right, bottom left, and bottom right panels 
show the distribution evolves after the BBHs have inspiralcd 
to r = 1000, 100 and 10M respectively, also as in Fig. 



the recoil velocity, while anti-alignment (cos #12 — > —1) 
does the opposite. 

The magnitude of this spin alignment is greatly re- 
duced for smaller mass ratios as seen in Fig. [5] for the 
case q = 1/3. Although the clustering of all the points 
about cos 6*12 = cosA</> is again apparent, the distribu- 
tions of the red, green, and blue points remain similar 
all the way down to r/ = 10M as seen in the lower right 
panel. The weaker influence of the spin-orbit resonances 
for g = l/3 follows from the smaller value of ri oc k in 
Eq. (|1.4[) . and is similarly reflected by the smaller frac- 
tion of the (cos 0i, cos 82 ) plane occupied by the resonant 
curves in Fig. [2J 

We have provided histograms of cos #12 and A<j> in 
Fig. [7] to clarify the differences between Figs. [5] and 
[H We see that the distributions of cos 6*12 and Acj) 
are initially flat for both mass ratios, but evolve con- 
siderably for q = 9/11 while remaining nearly flat for 
q = 1/3 within the limits set by Poisson fluctuations. 
The open blue (red) curves in the left panels of Fig. [7] 
clearly show distributions peaked at cos #12 = 1, A<fi = 0° 
(cos 0i2 = — 1, A0 = ±180°). Such trends are barely no- 
ticeable in the right panels. We will explore the impli- 



FIG. 6: Distributions of (A<f>, cos #12 ) at different separations 
r for 1000 initially isotropic maximally spinning BBHs with 
a mass ratio q = 1/3. The different panels, points, and lines 
are the same as those given for q = 9/11 in Fig. [5] 



cations of these findings for the final spins in the next 
Section. 



V. FINAL SPIN DISTRIBUTIONS 

Several attempts have been made to predict the final 
dimensionless spin %/ of the black hole resulting from a 
BBH merger. Initial attempts focused on finding sim- 
ple phenomenological fitting formulae for the final spin 
resulting from non-spinning, unequal-mass BBH merger 
simulations (4l| - [43| . A group at the Albert Einstein Insti- 
tute (AEI) developed a fitting formula that provides the 
magnitude and direction of y f in terms of the initial spins 
X\ , X2 an d the mass ratio q |44M46j . They assumed that 
the final spin magnitude could be expressed as a poly- 
nomial in xii X2, and the symmetric mass ratio r\, then 
made some additional assumptions about the symmetries 
of this polynomial dependence and how energy and an- 
gular momentum are radiated to reduce the number of 
terms in their expression. The coefficients of the remain- 
ing terms were calibrated using numerical-relativity (NR) 
simulations of BBH mergers in which the initial spins 
were either aligned or anti-aligned with the orbital angu- 
lar momentum. We shall refer to this older AEI formula 
as "AEIo" . A more recent paper [3(J by members of this 
group uses newer NR simulations to recalibrate their co- 
efficients, and replaces earlier assumptions with the con- 
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FIG. 7: Histograms of cos 812 and A(f> for BBHs with initially 
isotropic spins. The two left panels are for the mass ratio q = 
9/11, while the two right panels are for q — 1/3. The two top 
panels give the distribution of cos #12 , while the two bottom 
panels give the distribution of A<f>. The black curves are for all 
1000 BBHs in the 10 x 10 x 10 grid discussed in the text, while 
the blue (red) curves correspond to the blue (red) points in 
Figs. IMA with initial values cos6»i > 0.4 (cos#i < -0.4). The 
horizontal dotted lines show the initially flat distributions, 
while the solid lines show the distributions at r — 10Af. 



jecture that the final spin points in the direction of the 
total angular momentum of the initial BBH at any sepa- 
ration. For consistency, this requires the further assump- 
tion that angular momentum is always radiated in the 
direction of the total angular momentum. We shall refer 
to this newer AEI formula as "AEIn" . An alternative fit- 
ting formula was proposed by a group at Florida Atlantic 
University (FAU) [47| . Following the procedure outlined 
in [H,|39|, the FAU group performed 10 equal-mass mis- 
aligned simulations to calibrate the coefficients of fitting 
formulae for the Cartesian components of Xf ■ They then 
made additional assumptions about the mass-ratio de- 
pendence of these formulae, and found good agreement 
between their predictions and independent NR simula- 
tions with mass ratios as small as q = 5/8. We shall refer 
to the formula of this group as "FAU" . The Rochester In- 
stitute of Technology (RIT) group proposed yet another 
fitting formula during the preparation of this paper j48j . 
This formula includes higher-order terms in the initial 
spins that may ultimately be needed to describe future 
high-accuracy NR simulations. However, current simula- 
tions are inadequate to calibrate all the terms appearing 
in the RIT formula, so we will not consider its predictions 



in this paper. 

Other groups have predicted final spins by extrapo- 
lating analytical test-particle calculations to finite mass 
ratios, rather than calibrating fitting formulae with NR 
simulations. Buonanno, Kidder, and Lehner (BKL) [49[ 
derived a formula for the final spin by assuming, as is 
true in the test-particle limit, that the angular momen- 
tum radiated during the inspiral stage of a BBH merger 
exceeds that radiated during the plunge and ringdown. 
Using this assumption, they equated the final spin with 
the total angular momentum J = Lisco + Si + S2, where 
Lisco is the orbital angular momentum at the innermost 
stable circular orbit (ISCO) of a test particle of mass rjM 
orbiting a black hole of mass M and dimensionless spin 
Xf equal to that of the final black hole. This counterin- 
tuitive but inspired choice correctly provides Xf ~ I Xi m 
the q — > limit and respects the symmetry of BBH merg- 
ers under exchange of the labels of the two black holes. 
Though derived only from test-particle calculations, the 
BKL formula is remarkably successful at predictin g fi- 
nal spins even for equal-mass BBH mergers. Kesden [50| 
slightly modified the BKL spin formula to account for the 
energy radiated during the inspiral stage of the merger. 
This change makes the formula accurate to linear order 
in q in the test-particle limit. It gcncrically increases the 
magnitude of the predicted dimensionless final spin by 
reducing the predicted final mass m/ below M in the 
denominator of the expression Xf = Sfjm 2 ^. This in- 
crease improves the agreement with NR simulations of 
non-spinning BBH mergers, but leads to somewhat larger 
final spins than the other formulae for mergers of max- 
imally spinning BBHs, such as those considered in this 
paper. The predictions of this formula are refcred to as 
"Kcs" in this paper. 

We now present the predictions of the spin formulae 
summarized above for various distributions of BBH spins 
that are allowed to inspiral from = lOOOAf to rf = 
10M. 



A. Spin Magnitudes 

In the top panel of Fig. we show the final spin mag- 
nitude Xf predicted by the AEIn formula for the evenly 
spaced 10 x 10 x 10 grid of maximally spinning BBHs with 
q = 9/11 described in Sec. |IVl The other spin formulae 
give very similar results; the mean and variance of the 
final spin distributions predicted by the other formulae 
for some of the initial distributions described below are 
provided in Table fl] As in Figs. [3][71 the black curves 
in Fig. [8] refer to all 1000 BBHs, the blue curves to the 
subset of 300 BBHs with the lowest values of 6q, and 
the red curves to the subset of 300 BBHs with the high- 
est values of 9\. The dotted curves give the final spin 
distribution predicted for the BBH spin configurations 
at their initial separation n — 1000A/, while the solid 
curves give the final spin distribution predicted when 
these same BBHs arc allowed to inspiral to rf = 10M 
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FIG. 8: Top panel: Histogram of the final spin \f predicted 
by the AEIn formula for 1000 BBHs with mass ratio q = 
9/11 and isotropically distributed spins at rj = 1000 Af. The 
blue curves show the subset of 300 BBHs with the lowest 
initial values of 9i , while the red curves show the subset of 300 
BBHs with the highest initial values of 0\. The solid curves 
show the predicted spins if the AEIn formula is applied at 
rf — 10M after the BBHs have inspiraled to this separation 
according to the equations of Sec. [H] The dotted curves show 
the predicted spins if the AEIn formula is applied to the initial 
distribution at n = lOOOAf. Bottom panel: Histograms of the 
predicted final spins for 6 sets of BBH mergers with q = 9/11, 
and flat distributions in cos 82 and A<f> at n = 1000A/. The 
red, orange, yellow, green, blue, and purple curves have 9i = 
170°, 160°, 150°, 30°, 20°, and 10° respectively. As in the top 
panel the final spins predicted by applying the AEIn formula 
at Vi = lOOOAf are shown by dotted curves, while allowing the 
BBHs to inspiral to r/ = 10M before applying the formula 
leads to the spins shown by the solid curves. 



according to the PN evolution described in Sec. HH The 
AEIn formula is unique in that it claims to accurately 
predict final spins at all separations; separations as large 
as r = 2 x 10 4 M were considered in [3(J|. The other fit- 
ting formulae were intended to apply at tnr — 10M, the 
starting point for the NR simulations with which their 
coefficients were calibrated. The BKL and Kcs formu- 
lae were designed for use at the ISCO. Although strictly 
speaking the formulae other than AEIn cannot be ap- 
plied to widely separated BBHs, one can imagine that 
the BBHs inspiral to rf = 10M without spin precession 
where these formulae are valid. It is in this sense that we 
consider the predictions of these other formulae when we 
claim in this Section to apply them to BBH spin config- 
uration at n = 1000 M. 



FIG. 9: Histograms of the final spins \S predicted by the 
AEIn formula for the same sets of BBHs presented in Fig. [8] 
but with the mass ratio q — 1/3 instead of q — 9/11. As in 
that figure, the predictions made at = 1000A4" are shown 
with dotted curves, those made at r/ = 10M are shown with 
solid curves. The black curves in the top panel show the full 
set of 1000 BBHs, while the blue (red) curves show the subset 
of 300 BBHs with the lowest (highest) initial values of 9\. In 
the lower panel, the red, orange, and yellow curves show BBHs 
with Xi initially anti-aligned with Lat (#1 = 170°, 160°, 150°). 
The green, blue, and purple curves show BBHs with Xi m i~ 
tially aligned with Lat (<?i = 30°, 20°, 10°). 



The dotted and solid black curves in the top panel of 
Fig. [5] are identical to within the Poisson noise of our 
limited number of BBH inspirals, confirming the finding 
of Refs. [H, [H H3 that isotropic distributions of BBH 
spins remain nearly isotropic as they inspiral. Even at 
ri = lOOOM, the blue (red) subset of spin configura- 
tions yields the largest (smallest) predicted final spins, 
because for these configurations the spin of the more mas- 
sive black hole is aligned (anti-aligned) with the orbital 
angular momentum. The spin-orbit resonances further 
enhance (reduce) the final spins predicted for these sub- 
sets by aligning (anti-aligning) the BBH spins with each 
other during the inspiral for small (large) initial values 
of 6\ . As a result, the solid blue (red) distribution at 
rf = 10AI has a larger (smaller) mean final spin than 
the initial dotted distribution at r, = 1000 M. This can 
be seen in the displacement of predicted final spins for 
the colored subsets away from \f — 0-75 towards larger 
and smaller values. 

To clarify the magnitude of this effect, we have per- 
formed 6 additional sets of BBH inspirals, each of which 
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consists of a fixed value of 6\ and a 30 x 30 grid evenly 
spaced in cos 82 and Ac/). Three of these sets have the 
spin of the more massive black hole nearly aligned with 
the orbital angular momentum {6\ = 10°, 20°, 30°), while 
the other 3 sets have Xi nearly anti-aligned with L^v 
(01 = 150°, 160°, 170°). The choice of aligned distribu- 



tions was partly motivated by the finding of Ref. [17 1 
that accretion torques will align BBH spins to within 10° 
(30°) of the orbital angular momentum for a cold (hot) 
disk. The predicted final spins for these distributions, 
both at ri = 1000 M and rf = 10M, are shown in the 
bottom panel of Fig. [5] The final spins for the initially 
aligned (6q < 30°) BBH distributions are significantly 
larger when predicted at rj_ = 10M than at = 1000 M . 
undermining the claim of [30J that the AEIn formula can 
accurately predict final spins at large separations with- 
out the need for PN evolutions. The predicted final spins 
for the initially anti-aligned (81 > 150°) BBH distribu- 
tions conversely shift to lower values as the predictions 
are made later in the inspiral. We provide the mean and 
standard deviation of the final spins predicted for these 
6 new sets of partially aligned BBH distributions for all 
5 formulae in Table HI 

To explore the dependence of these effects on the mass 
ratio, we have provided histograms of the predicted final 
spins for these same BBH spin distributions with q = 
1/3 in Fig. [S] The discrete peaks at low values of Xf m 
the histograms in the top panel are an artifact of the 10 
discrete values of cos 81 in our 10 x 10 x 10 grid. Each peak 
contains 100 points with the same initial value of 8\. The 
decrease in the width of each peak as the BBHs inspiral 
from n = 1000M to rf = 10M is a consequence of the 
anti-alignment of the BBH spins for large 8\ , but the gaps 
between the peaks would be filled in if we used a finer 
grid. The shifts in the mean values of the peaks should 
be robust with respect to the grid spacing. These shifts 
for the initially aligned BBH distributions are provided 
in Table U for all 5 formulae for q = 1/3, as well as for 
the intermediate mass ratio q = 2/3. 



B. Spin Directions 

Before providing quantitative results, we need to clar- 
ify what is meant by the direction of the spin of the final 
black hole. In what reference frame is this direction de- 
fined? Most of the fitting formulae calibrated with NR 
simulations attempt to predict the angle 



f 



(5.1) 



between the BBH orbital angular momentum Ljv at the 
separation r j = tnr where the NR simulations were per- 
formed and the final spin Xf predicted from the BBH 
spin configuration at this same separation. The analyti- 
cal predictions of BKL and Kes were designed to apply to 
BBH spin configurations at rf = nsco- If one assumed 
that neither the orbital angular momentum nor the BBH 



spins (upon which the prediction Xf( r f) depends) pro- 
cessed during the inspiral, one could insert these quanti- 
ties at any separation into the right-hand side of Eq. (|5.1[) 
to predict The angle is physically interesting be- 
cause it quantifies the post-merger alignment between x / 
and the inner edge of the accretion disk if one assumes 
that torques have aligned the circumbinary disk with Ljv- 
However, one might also be interested in the alignment 
between Xf an d a feature like the galactic disk that is 
assumed to be aligned with Ljy at some larger scale r^. 
In that case, one would need to compute the angle 



$i = arccos[Ljv(ri) • Xf(n)] 



(5.2) 



between Ljv at this larger separation and the final spin 
Xfi r i) predicted from the BBH spins at this same sepa- 
ration. 

The proper way to predict % f from the BBH spins at 
r,; would be to use PN equations like those specified in 
Sec. [H]to propagate those spins and Ln to down to r/, 
then insert them into the fitting formula of one's choice. 
The AEIn formula is based on the conjecture that Xf 
points in the direction of the total angular momentum 
J at any separation, since angular momentum is always 
radiated parallel to J, thus preserving its direction. This 
conjecture is plausible because at large separations, the 
precession time t p is much shorter than the inspiral time 
tow- If the vectors associated with the BBHs precess 
rapidly enough, all components except those parallel to J 
(which varies on the longer timcscale £gw) will average to 
zero. The AEIn conjecture is very useful because it allows 
to be computed without solving any PN equations. 
However, the approximation t p <C icw upon which it 
depends breaks down at small separations. This may 
lead to incomplete cancellation of the angular momentum 
radiated perpendicular to J. 

We test this possibility by calculating 



8 j = arccos[J(ri) • J(ry)] , 



(5.3) 



the angle between the total angular momentum at r,; = 
1000M and that after the BBHs have inspiraled to rf = 
10M. If the direction of J really was preserved during the 
inspiral, 8j would vanish. We present histograms of 8 j 
for mass ratio q = 9/11 in Fig. 1101 The upper panel shows 
the 10 x 10 x 10 grid of BBH spin configurations evenly 
spaced in (cos 8\ , cos 82 , At/>) that we have discussed pre- 
viously. The direction of J changes by 8j < 2° during 
most of the inspirals, though a tail extends to larger val- 
ues for large initial values of 8\. This tail can be seen 
more clearly in the bottom panel for the BBHs with Xi 
initially anti-aligned with (6>i > 150°). We agree 
with [30] that these large changes in the direction of J 
are likely a consequence of the transitional precession first 
identified in Ref. [5l| . This transitional precession occurs 
to an even greater extent for smaller mass ratios, as can 
be seen in Fig. [IT] for q = 1/3. As in the upper panel of 
Fig. [5J discrete peaks resulting from the grid spacing in 
cos 8i can be seen in the left panel of Fig. [TTJ The mid- 
dle panel shows that the direction of J remains nearly 
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FIG. 10: Top panel: Histogram of the angle 9j (in degrees) 
between the total angular momentum J at n = 1000M and 
that at 77 = 10M for our set of 1000 BBHs with q = 9/11 
and initially isotropic spins. As in previous figures, the blue 
(red) curve shows the subset of 300 BBHs with the lowest 
(highest) initial values of 81 . Bottom panel: Histograms of 
8 j for the 6 sets of 900 BBH mergers with flat distributions 
in cos 82 and Acj> at = 1000M. The red, orange, yellow, 
green, blue, and purple curves show BBHs that have 8\ — 
170°, 160°, 150°, 30°, 20°, and 10° respectively at this initial 
separation. 



constant (9j < 0.5°) when Xi m closely aligned with 
Ljv (#i < 30°). However, the right panel shows that 
the assumption of constant J fails badly for the BBHs 
with 6*i > 150°, that comprise ~ 7% of isotropically dis- 
tributed BBH mergers. The mass ratio q = 1/3 is not 
extreme compared to the majority of astrophysical merg- 
ers, so caution should be taken when assuming that Xf 
points in the direction of J such as in Eq. (|5.2p . 

What about the less ambitious predictions of from 
BBH spins at 77 = 10M, assuming that NR simulations 
correctly describe spin precession from this separation 
until merger? Spin-orbit resonances have significant im- 
plications for these predictions as well. We show pre- 
dictions of i?/ by the AEIn formula for a mass ratio of 
q = 9/11 in Fig. 1121 The other formulae predict very sim- 
ilar results. As in Figs. E] and O the dotted curves show 
predictions assuming that the initial BBH spin distribu- 
tion is preserved down to rj = 1QM . The solid curves 
include spin precession from = lOOOAf to 77 = 10M 
according to the PN equations of Sec. HU The differ- 
ence between the dotted and solid black curves in the 
top panel is below the Poisson fluctuations, another con- 



FIG. 11: Left panel: Histogram of the angle 8j (in degrees) 
between the total angular momentum J at n — lOOOAf and 
that at 77 = 10M for our set of 1000 initially isotropically 
spinning BBHs with q — 1/3. As in previous figures, the blue 
(red) curve shows the subset of 300 BBHs with the lowest 
(highest) initial values of 0\. Middle panel: Histograms of 8j 
for the 3 sets of 900 BBH mergers initially with 81 = 10° (pur- 
ple), 20° (blue), and 30° (green). Right panel: Histograms of 
8 j for the 3 sets of 900 BBH mergers initially with 9i = 150° 
(yellow), 160° (orange), and 170° (red). 

sequence of the finding of Refs. [23|, H(|[22} that isotropi- 
cally oriented BBH spins remain nearly isotropic as they 
inspiral. Careful examination of the upper panel reveals 
that spin precession has shifted the BBHs with Xi m i~ 
tially aligned with (blue distribution) to larger -&f, 
while the anti-aligned BBHs have conversely shifted to 
smaller fif. 

This trend is much more pronounced in the middle and 
bottom panels of Fig. [T3] Spin precession actually results 
in the initially aligned BBHs {6\ < 30°) having larger 
values of at 77 = 10M than the anti-aligned BBHs 
(9i > 150°), a reversal of what would be predicted from 
the initial spin distributions shown by the dotted curves. 
The spin-orbit resonances explain this highly counterin- 
tuitive result. The BBHs initially with 6\ < 30° are 
influenced by the Acf> — 0° resonances which align the 
BBH spins with each other and anti- align S = Si + S2 
with Lat. Both effects lead to larger predicted values of 
■df. Conversely, the BBHs initially with 6\ > 150° are 
influenced by the Acj> = 180° resonances, which greatly 
decrease the magnitude of S and align it with Ljv- This 
explains the reduced values of 1?/ for these BBHs seen 
in the bottom panel of Fig. [T2] This same effect can 
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FIG. 12: Top panel: Histogram of the angle (in degrees) 
between the orbital angular momentum Lat at 77 = 10M and 
the final spin Xf predicted by the AEIn formula from the 
BBH spins at that separation. The BBHs have a mass ratio 
q = 9/11. As in previous figures, the black curves show 1000 
mergers with initially isotropic BBH spins, while the blue 
(red) curves show the subset of 300 BBHs with the lowest 
(highest) initial values of Q\. The dotted curves show predic- 
tions in the absence of spin precession, while the solid curves 
show how these predictions change when the BBH spins pre- 
cess from = 1000 M to r/ = 10M according to the PN 
equations of Sec. [II] Middle panel: Histograms of for the 
3 sets of 900 BBH mergers with Xi initially aligned with Lat 
[0i = 10° (purple), 20° (blue), 30° (green)]. Bottom panel: 
Histograms of i?/ for the 3 sets of 900 BBH mergers with 
Xi initially anti-aligned with Ljv [61 = 150° (yellow), 160° 
(orange), 170° (red)]. 

be seen for a mass ratio of q = 1/3 in Fig. Q2] albeit 
with less significance owing to the weaker resonances at 
this smaller mass ratio. Figs. [T^] and Q2] again illustrate 
the importance of accounting for spin precession between 
r.; = 1000M and 77 = 10M when attempting to predict 
final spins. 

VI. SPIN PRECESSION UNCERTAINTY 

So far, we focused on how spin precession between 
and 77 alters the expected distribution of final spins. In 
this Section, we show that spin precession introduces a 
fundamental uncertainty in predicting the final spin. An 
uncertainty A7' in the BBH separation leads to an uncer- 
tainty Atcw m the time until merger. If this uncertainty 
is comparable to the precession time t pi the phase of the 



FIG. 13: Histograms of the angle predicted by the AEIn 
formula for the same sets of BBHs shown in Fig. [12] but with 
a mass ratio q = 1/3. As in that figure, the top panel shows 
BBHs with initially isotropic spins with the blue (red) curves 
indicating those BBHs with the lowest (highest) initial values 
of 81 . Dotted curves show predictions without spin precession, 
while the solid curves show how these predictions change if the 
BBHs spins precess from n — 1000 M to 77 = 10M according 
to the PN equations of Sec. [H] The bottom panel shows 
distributions with flat initial distributions of 0036*2 and A0, 
but with 0i now initially set to 170°, 160°, 150°, 30°, 20°, and 
10° respectively for the red, orange, yellow, green, blue, and 
purple curves. 

spin precession at which the merger occurs will be un- 
certain as well. This new uncertainty is independent of 
and may exceed that associated with the NR simulations 
themselves. Readers only interested in astrophysical dis- 
tributions of final spins may wish to proceed to the dis- 
cussion in Sec. IVIII 

It is often useful to define the final spin direction rel- 
ative to the orbital angular momentum Ljv at different 
separations. We therefore generalize the angles defined 
in Eqs. (|5.1[) and ()5.2[) to the separation-dependent quan- 
tities 

i?/(r) = arccos[Ljv(r) -X/(r)] , (6.1) 
#i(r) = arccos[Liv(ri) • Xf(r)] . (6.2) 

Note that these quantities reduce to the previously de- 
fined angles in the appropriate limit: $/(/■/) = 
^i(fi) = These definitions address two ambiguities; 
(i) the choice of the reference orbital angular momentum 
and (ii) the separation at which a given fitting formula 
is evaluated. 
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FIG. 14: Predicted Xf (upper panels) and ~&f{r) (lower panels) obtained from the AEIn (left) and the Kesden (right) formula. 
The initial parameters of the binary are q — 9/11, xi — X2 — 1, #i = 120°, 62 = 60° and A(f> — 288°. The different curves 
correspond to initial frequencies Mwi = 3.16 x lCT 5 (solid), 3.17 x 1CT 5 (long-dashed), 3.18 x 1(T 5 (dashed) and 3.19 x 1(T 5 
(dotted). The envelope determined for Muj = 3.16 x 1CP 5 is displayed by thin solid curves. The upper horizontal axis gives 
the binary separation in units of M; the lower horizontal axis gives the corresponding orbital frequency Afcj. 



Before we discuss the uncertainties in determining 
these angles and the final spin magnitude, we illustrate 
the evolution of these quantities during the PN inspi- 
ral for a few characteristic examples. In Fig. [TJ] we 
display the final spin magnitude x/( r ) an d the angle 
t?/(r) as predicted by the AEIn and the Kesden for- 
mulae for a binary with mass ratio q — 9/11, extremal 
spins, and initial spin orientation specified by the an- 
gles 0i = 120°, 2 = 60°, A(j) = 288°. The behavior 
of the AEIo, FAU and BKL formulae is quite similar to 
the Kesden formula. The different curves in each panel 
correspond to slightly different initial frequencies or sep- 
arations, Mui = 3.16 x 10~ 5 , 3.17 x 1CT 5 , 3.18 x 10~ 5 
and 3.19 x 10 . The spin precession gcnerically mani- 
fests itself in the oscillatory character of the curves; these 
oscillations would be absent for the resonant configura- 
tions described in Sec. IIII1 The thin solid lines represent 
envelope functions obtained by fitting fourth-order poly- 



nomials to the maxima and minima, respectively, of the 
evolutions starting with Mu>i = 3.16 x 10~ 5 . Note that 
these fits contain no information on the results obtained 
by using different values of Muji, and yet they still pro- 
vide excellent envelopes in all cases. 

This figure illustrates two ambiguities in predicting x / : 
(i) the initial frequency uii at which the BBH parameters 
are specified, and (ii) the final separation 77 at which the 
given formula for Xf should be applied. Uncertainty in 
the separation at which the binary decouples from exter- 
nal interactions could lead to ambiguity in oji in theoret- 
ical studies, while uncertainty in the observed distance, 
projected separation, or line-of-sight velocity could lead 
to uncertainty in cjj for models of particular systems. 
Gauge-dependent definitions of rj could lead to uncer- 
tainty in the separation at which fitting formulae should 
be applied. Our task in evaluating the resulting uncer- 
tainties for the fitting formulae AEIn, AEIo, FAU, BKL 
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FIG. 15: Predicted spin direction i?;( r ) from the AEIn (left panels) and the Kes (right panels) formula. The upper panel shows 
the evolution of a binary starting with initial parameters g = 9/11, xi — X2 = 1, 61 = 120°, 82 = 60° and A(f> = 288°, as in 
Fig. [14] For comparison, in the bottom panels we consider a mass ratio q = 1/3 and initial spin parameters \i — X2 — 1> 
0i = 154°, 8 2 = 124° and A0 = 216°. The initial frequency is Muji = 3.16 x 10~ 5 (solid), 3.17 x 10" 5 (long-dashed), 
3.18 x 10" 5 (dashed) and 3.19 x 10" 5 (dotted curve). The upper horizontal axis gives the binary separation in units of M; the 
lower horizontal axis gives the corresponding orbital frequency Moj. 



and Kes introduced in Sec. [V] is somewhat simplified be- 
cause both ambiguities are rooted in the rapid variations 
of the phase and in the resulting oscillations in the final 
quantities. These precession-induced oscillations are a 
clear manifestation of the hierarchy of time scales intro- 
duced in Eq. ([53]) : t p < t GW . 

In the upper panels of Fig. 1151 we show the angle 
for the same binary configuration illustrated in Fig. 1141 
In the lower panel of Fig. [TS] we consider instead, for 
comparison, a system with lower mass ratio (7 = 1/3 and 
initial spin orientation 6 X = 154°, 8 2 = 124°, A</> = 216°. 
As before, different curves correspond to different initial 
frequencies. The predicted spin direction as described by 
$j(r) shows little variation with w,-. On the other hand, 
the figure demonstrates a strong dependence of on 
the separation r at which we apply the fitting formulae. 

In the remainder of this Section, we discuss the uncer- 



tainties caused by the rapid spin precession of the follow- 
ing quantities: 

1) X/( r /) : The magnitude of the final spin as pre- 
dicted by applying a given fitting formula at small 
binary separation rf = 10M, i.e. shortly before 
merger. 

2) $f(rf): The angle between the orbital angular mo- 
mentum at rf = 10M and the final spin (as pre- 
dicted using the binary parameters at rj = 10 M). 

3) z?i(r/): The angle between the orbital angular mo- 
mentum of the binary at large separation and the 
final spin (as predicted using the binary parameters 
at 77 = 10M). 

4) $i(r"j): The angle between the orbital angular mo- 
mentum of the binary at large separation and the 
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FIG. 16: Uncertainties in the final spin magnitude Xf( r f) (left) and direction (right) for extremal BBHs with mass 

ratio q = 2/3. Solid black histograms were obtained by starting the evolutions at Muji = 3.16 x 10~ 5 and using the envelope 
method of Fig. [14] Dashed red histograms were obtained by considering the maximum variation in the final quantities as we 
let Mui vary from 3.16 x 1CT 5 to 3.22 x 1(T 5 in steps of 0.005 x 10~ 5 . 



final spin predicted using the binary parameters 
at this same large separation. We investigate the 
claim that the AEIn formula, unlike the others, can 
determine this angle without evolving the BBH pa- 
rameters down to Tj. 

These quantities are important for modeling the as- 
sembly of supermassive black holes in the context of 
cosmological structure formation (see e.g. 0, I521456T ] ) . 
They are also relevant for electromagnetic counterparts 
of gravitational- wave sources [57j , especially when the in- 
voked mechanism producing the counterparts depends on 
the recoil velocity of the remnant black hole [58l460| |. 

We determine the precession-induced uncertainties as 
follows. Individual evolutions, such as those considered 
in Fig. [TJJ suggest that the width of the envelopes or the 
dispersion induced by varying the initial frequency pro- 
vide very similar estimates for the uncertainty in X/( r /) 
and $/(r/). We have verified this conjecture by evolving 
the evenly spaced 10 x 10 x 10 grid of initially isotropic, 
maximally spinning BBH configurations introduced in 
Sec. I1VI for mass ratio q = 2/3 and several slightly dif- 
ferent initial frequencies. When we estimate uncertain- 
ties by varying Moj 1 from 3.16 x 10~ 5 to 3.22 x 10~ 5 
in steps of 0.005 x 10 -5 we obtain the red dashed his- 
tograms in Fig. [T5] These histograms are in good agree- 
ment with the black solid histograms, where the uncer- 
tainty was estimated from the width of the envelopes. In 
order to reduce computational cost, in the remainder of 
this Section we determine the uncertainties A%/(r/) and 
Ai)f(rf) by evolving an ensemble of binaries from a sin- 
gle initial frequency (Muii = 3.16 x 10 -5 ) and using the 
envelope method. 

Fig. [T5] shows that the envelope method does not ade- 
quately describe the uncertainty in "dn (r) . Why docs this 
angle behave so differently from $/(r) as illustrated in 



Fig. 1141 .'' The direction of Ljv(rj) is fixed, while accord- 
ing to the AEIn formula Xf( r ) points in the direction 
of J(r). The total angular momentum J only varies on 
the radiation timescale 4gW) so according to Eq. (|6.2[) 
the AEIn prediction of should only vary on this 

slower timescale as well. The left panels of Fig. [T51 at 
least at small orbital frequencies Muj where t p <C tow, 
indeed lack the high-frequency oscillations characteristic 
of spin precession. In contrast, the Kesden predictions for 
$j(r) shown in the right panels of Fig. [T5] are varying in a 
more complicated way than the predicted values of $f(r). 
Changes in the angle $i(r) between the fixed L/v(?*i) and 
varying Xf( r ) reflect the full complexity of spin preces- 
sion for misaligned, unequal-mass BBHs. The simpler 
variation in #/(r) occurs because both Ljv(r) and Xf( r ) 
are jointly precessing about 3(r), albeit on the same short 
timescale t p . 

Since the envelope method fails for $i(r), we some- 
what arbitrarily define the uncertainty A #3(77) as the 
maximum deviation of ??j(r) from "di{rf = 10M) in the 
window Tf < r < 2r/. This window covers approxi- 
mately the range of initial separations within the reach 
of present and near-future numerical relativity simula- 
tions, while smaller separations must be excluded due to 
the breakdown of the PN expansion. The formulae other 
than AEIn do not claim to predict x / from the BBH pa- 
rameters at large separations. To apply these formulae 
correctly, one must evolve the BBH parameters inwards 
to rj according to PN equations such as those in Sec. |TT] 
before applying the formulae. This evolution requires sig- 
nificant additional effort, but if performed properly would 
only increase Ai?,(?"j) above Ai?j(r/) by the uncertainty 
in the PN equations themselves. The uncertainty com- 
ing from PN evolutions could be quantified by comparing 
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different PN orders and pushing the calculation of spin 
contributions to higher order; such an analysis is beyond 
the scope of this paper. The AEIn formula is special in 
that it predicts fl^ri) without this additional PN evo- 
lution. Since AEIn claims that both Ljv(ri) and Xf( r ) 
are independent of r, the uncertainty A$j(r,) for this 
formula is the maximum deviation from $i(r/ = 10M) 
over the entire interval < r < Since the orbital 
angular momentum L jv increasingly dominates over spin 
contributions in the sum J = L^v + Si + S2 at large sep- 
arations, Lat has little opportunity to precess at large 
separations and the uncertainty i?i(r,) asymptotes to a 
constant value in this limit. 

We have evolved the uniform 10x10x10 grid of max- 
imally spinning binaries introduced in Sec. IIVI for three 
different mass ratios: q = 9/11, q = 2/3 and q = 1/3. 
The average uncertainties (plus or minus their associated 
standard deviations) are summarized in Table [TTJ 

Errors in the final spin magnitudes due to the rapid 
spin precession are in the range Axf < 0.03 for all mass 
ratios. The FAU formula performs exceptionally well for 
nearly equal masses, although it deteriorates to the level 
of the other predictions for q = 1/3. We suspect that this 
is because several of the higher-order terms in r\ in the 
FAU formula are symmetric in the dimcnsionless spins 
Xij X21 while physically one would expect the spin of the 
more massive black hole to be more important in the 
limit q — > 0. Overall however, all formulae are able to 
predict the spin magnitude with rather good accuracy. 

The uncertainty (r/) in the angle between the fi- 
nal spin and the orbital angular momentum shortly be- 
fore merger is typically in the range of a few to 20 de- 
grees. Investigation of the angular dependence of the 
spin uncertainties shows that the AEIn formula tends to 
behave better for initially aligned spins (small 9\ and 62 ) 
and worse for anti-aligned cases. This is likely a conse- 
quence of anti-aligned binaries being closer to the limit 
L(r) m — S(r) where transitional precessio n 1511 occurs, 
violating assumptions (iii) and (iv) of Ref. [30( . 

All formulae are able to predict the angle $i(r/) be- 
tween the initial orbital angular momentum and the fi- 
nal spin with decent accuracy. The AEIn predictions are 
overall more accurate, but investigation of the angular 
dependence reveals that this accuracy deteriorates (as 
expected) when q = 1 /3 and the spin of the larger black 
hole is nearly anti-aligned. In this limit the uncertain- 
ties increase up to ~ 20°. This is again a consequence of 
those configurations approaching the transitional preces- 
sion regime, where L(r) ~ — S(r). 

The AEIn prediction is unique in that it claims to pre- 
dict ^j(rj) using the binary parameters at large separa- 
tion without PN evolution. Our findings confirm (quite 
remarkably) that the majority of binaries in an initially 
isotropic ensemble result in a final spin which is nearly 
aligned with the orbital angular momentum at large bi- 
nary separation. The values of 9j shown in Figs. HU1 and 
[TT] suggest that this would not be the case for BBHs ini- 
tially anti-aligned with Ljv- The accuracy of the AEIn 



predictions also decreases for unequal masses (as ex- 
pected and verified by our results for q = 1/3). More 
extreme mass ratios are expected to play a significant 
and possibly dominant role in the coalescence of SBH bi- 
naries |6l| - l63l ] , so it will be crucial to test the robustness 
of the Barausse-Rezzolla predictions for q = 1/10 and be- 
yond. Accurate PN evolutions are more difficult in this 
regime, and we plan to investigate more extreme mass 
ratios in the future. 



VII. DISCUSSION 

In this paper, we examined how precession affects the 
distribution of spin orientations as BBHs inspiral from an 
initial separations r, ~ 1000A/ where gravitational radi- 
ation begins to dominate the dynamics, all the way down 
to separations r/ — 10-M where numerical- relativity sim- 
ulations typically begin. 

We confirmed previous findings that isotropic spin dis- 
tributions at r,; ~ 1000 M remain isotropic at r/ ~ 10M 
[23l I2H l27l ]. However, torques exerted by circumbinary 
disks may partially align BBH spins with the orbital 
angular momentum at separations r > r% before gravi- 
tational radiation drives the inspiral [23|]. Recent sim- 
ulations suggest that the residual misalignment of the 
BBH spins with their accretion disk could typically be 
~ 10°(30°) for cold (hot) accretion disks, respectively 
Partially motivated by these findings, we carried 
out a more careful analysis of spin distributions that are 
partially aligned with the orbital angular momentum at 
r = i\. We found that spin precession efficiently aligns 
the BBH spins with each other when the spin of the more 
massive black hole is initially partially aligned with the 
orbital angular momentum, increasing the final spin. We 
found the opposite trend when the spin of the more mas- 
sive black hole is initially anti-aligned with the orbital 
angular momentum. Long evolutions are necessary to 
capture the full magnitude of the spin alignment. This 
could explain why these trends were not observed in the 
PN evolutions by Lousto et al. [27| . which began at a 
fiducial binary separation r = 50M. 

Some models of BBH evolution (see e.g. [Il|,[62j]) sug- 
gest that SBH mergers might have comparable mass ra- 
tios (q < 1) at high redshift and more extreme mass ra- 
tios at low redshift. Since spin alignment is stronger for 
comparable-mass binaries, more alignment might be ex- 
pected in SBH binaries at high redshifts. Observational 
arguments (see e.g. [641 ]) and magnetohydrodynamic sim- 
ulations of accretion disks [65j provide some evidence 
that black hole spins are related to the radio loudness 
of quasars. If so, the inefficient alignment (and conse- 
quently smaller spins) produced by unequal-mass merg- 
ers at low redshift would at least be consistent with recent 
observational claims that the mean radiative efficiency of 
quasars decreases at low redshift [6(1 [67| ■ Stellar- mass 
black hole binaries should also have comparable mass ra- 
tios, so significant spin alignment could occur in such 
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systems as well. 

We also pointed out that predictions of the final spin 
Xj usually suffer from two sources of uncertainty: (i) 
the uncertainty in the initial frequency Wi at which the 
BBH parameters are specified, and (ii) the uncertainty 
in the final separation 77 at which the given formula for 
Xj should be applied. Both ambiguities are rooted in the 
rapid preccssional modulation of the orbital parameters, 
which in turn results from the precessional timescale t p 
being much shorter than the radiation timescale £gw- 
Spin precession induces an intrinsic inaccuracy A%f < 
0.03 in the dimensionless spin magnitude and A-df < 20° 
in the final spin direction. 

The spin-orbit resonances studied in this paper should 
have significant effects on the distribution of gravitational 
recoil velocities resulting from BBH mergers, because the 
maximum recoil velocity has a strong dependence on spin 
alignment [l5l - fl7j . We plan to extend this study to inves- 
tigate the predictions of different formulae for the recoil 
velocities that have been proposed in the literature. 
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TABLE I: Mean and standard deviation of the final spin magnitudes predicted for different sets of maximally spinning BBH 
mergers. The first column lists the formulae used to predict the final spins, as described in Sec. [V] The second column gives 
the mass ratio q. Each set of BBHs begins at n — lOOOAf with the indicated value of 0i and flat distributions of cos 02 and A<f). 
The third, fourth, and fifth columns show the mean and deviation expected if the BBH spins do not precess, thus maintaining 
their initial distributions at n = lOOOAf until merger. The sixth, seventh, and eighth columns assume that the spins precess 
according to the PN equations of Sec. [II] as they inspiral to rf — 10M, at which separation we apply the spin formulae. 









n = lOOOAf 






r f = 10M 




model 


q 


0i = 10° 


0i = 20° 


6>i = 30° 


0i = 10° 


0i = 20° 


0i = 30° 


AEIn 


9/11 


0.867 ±0.064 


0.863 ±0.065 


0.857 ± 0.066 


0.914 ±0.034 


0.905 ± 0.036 


0.892 ± 0.038 


AEIo 


9/11 


0.866 ± 0.063 


0.863 ±0.064 


0.856 ± 0.065 


0.912 ±0.034 


0.904 ± 0.036 


0.891 ± 0.038 


FAU 


9/11 


0.873 ± 0.059 


0.868 ±0.060 


0.861 ± 0.061 


0.909 ± 0.035 


0.901 ± 0.037 


0.888 ± 0.039 


BKL 


9/11 


0.862 ± 0.067 


0.858 ±0.068 


0.851 ±0.070 


0.905 ± 0.037 


0.898 ± 0.039 


0.884 ± 0.042 


Res 


a / 1 1 
9/ 11 


a aai _i_ a ato 
0.901 ± 0.0 1 2 


0.89b ± 0.0/3 


a oon _l a n7c 
0.889 ± 0.0/0 


A AKA _L A AOO 

0.950 ± 0.038 


A A/1 1 _L A A/1 1 

0.941 ± 0.041 


A H07 _L A A/1 A 

0.92 / ± 0.044 


AEIn 


2/3 


0.886 ± 0.052 


0.882 ± 0.053 


0.875 ± 0.054 


0.922 ± 0.030 


0.914 ± 0.031 


0.900 ± 0.034 


AEIo 


2/3 


0.886 ± 0.052 


0.882 ± 0.052 


0.876 ± 0.054 


0.922 ± 0.030 


0.914 ± 0.031 


0.900 ± 0.034 


FAU 


2/3 


0.901 ± 0.043 


0.895 ± 0.044 


0.886 ± 0.046 


0.924 ± 0.029 


0.915 ± 0.030 


0.901 ± 0.031 


BKL 


2/3 


0.882 ± 0.052 


0.878 ± 0.053 


0.870 ± 0.054 


0.914 ± 0.031 


0.906 ± 0.032 


0.893 ± 0.035 


Kes 


2/3 


0.921 ± 0.056 


a m 7 i a ac^t 

0.917 ± 0.057 


a f\f\f\ 1 a ntrn 

0.909 ± 0.059 


A A C O 1 AAOI 

0.958 ± 0.031 


AA/1A J A AO /I 

0.949 ± 0.034 


A HOC L A AOT 

0.935 ± 0.037 


AEIn 


1/3 


0.950 ± 0.025 


0.946 ± 0.025 


0.938 ± 0.026 


0.957 ± 0.023 


0.951 ± 0.023 


0.941 ± 0.022 


AEIo 


1/3 


0.958 ± 0.025 


0.953 ±0.026 


0.944 ± 0.026 


0.964 ± 0.023 


0.958 ± 0.023 


0.947 ± 0.022 


FAU 


1/3 


0.972 ± 0.013 


0.964 ± 0.014 


0.951 ± 0.016 


0.975 ± 0.012 


0.966 ± 0.012 


0.953 ± 0.011 


BKL 


1/3 


0.931 ± 0.020 


0.927 ± 0.020 


0.921 ± 0.021 


0.936 ± 0.018 


0.931 ± 0.018 


0.923 ± 0.018 


Kes 


1/3 


0.968 ± 0.021 


a A£?tr i a noo 

0.965 ± 0.022 


A ACO 1 A AOO 

0.958 ± 0.023 


A AT A 1 A A 1 A 

0.974 ± 0.019 


A ATA 1 A A 1 A 

0.970 ± 0.019 


A A£?0 1 A AAA 

0.962 ± 0.020 


model 


q 


0i = 150° 


61 = 160° 


0i = 170° 


0i = 150° 


0i = 160° 


0i = 170° 


AEIn 


9/11 


0.551 ±0.080 


0.527 ±0.080 


0.511 ±0.079 


0.535 ± 0.072 


0.510 ±0.076 


0.493 ± 0.080 


AEIo 


9/11 


0.551 ±0.080 


0.527 ±0.080 


0.512 ±0.079 


0.535 ± 0.072 


0.510 ±0.077 


0.493 ± 0.080 


FAU 


9/11 


0.542 ± 0.076 


0.520 ±0.076 


0.506 ± 0.076 


0.530 ± 0.070 


0.507 ± 0.074 


0.492 ± 0.076 


BKL 


9/11 


514 ± 088 


488 ± 087 


471 + 086 


496 ± 078 


468 ± OS 1 } 


449 ± 087 


Kes 


9/11 


0.531 ±0.091 


0.504 ±0.090 


0.486 ± 0.089 


0.512 ± 0.081 


0.483 ± 0.087 


0.463 ± 0.091 


AEIn 


2/3 


0.500 ± 0.067 


0.467 ±0.066 


0.445 ± 0.065 


0.490 ± 0.057 


0.456 ± 0.062 


0.432 ± 0.065 


AEIo 


2/3 


0.499 ± 0.067 


0.466 ± 0.067 


0.444 ± 0.066 


0.489 ± 0.057 


0.455 ± 0.062 


0.432 ± 0.065 


FAU 


2/3 


0.490 ± 0.060 


0.460 ± 0.060 


0.441 ± 0.059 


0.483 ± 0.053 


0.452 ± 0.056 


0.432 ± 0.059 


BKL 


2/3 


0.465 ± 0.072 


0.430 ±0.071 


0.405 ± 0.070 


0.454 ± 0.060 


0.416 ±0.066 


0.390 ± 0.070 


Kes 


2/3 


0.480 ± 0.075 


0.442 ±0.071 


0.417 ±0.072 


0.468 ± 0.063 


0.428 ± 0.068 


0.401 ± 0.072 


AEIn 


1/3 


0.324 ± 0.034 


0.233 ±0.034 


0.151 ±0.032 


0.323 ±0.013 


0.231 ±0.015 


0.145 ±0.021 


AEIo 


1/3 


0.321 ± 0.034 


0.230 ±0.034 


0.146 ±0.032 


0.319 ±0.012 


0.227 ±0.015 


0.140 ±0.021 


FAU 


1/3 


0.301 ± 0.026 


0.222 ±0.026 


0.154 ±0.024 


0.300 ± 0.016 


0.220 ±0.018 


0.151 ±0.020 


BKL 


1/3 


0.315 ± 0.034 


0.222 ±0.034 


0.136 ±0.032 


0.313 ±0.011 


0.219 ±0.013 


0.130 ±0.019 


Kes 


1/3 


0.322 ± 0.035 


0.227 ±0.035 


0.139 ±0.033 


0.320 ±0.012 


0.224 ±0.014 


0.133 ±0.019 



TABLE II: Uncertainty distributions in \f an d in the various angles describing the final spin directions, as predicted by the 
formulae listed in Sec.[V] The uncertainties and their standard deviations are obtained by evolving uniform 10 x 10 x 10 grids 
of maximally spinning BBHs with mass ratio q — 9/11, 2/3 and 1/3, respectively. 



model 


q 


A X f(r = 10M) 


Ad f (r = 10 M) 


A#i(r = 10M) 


Ai?,(r = lOOOAf) 


AEIn 


9/11 


0.0159 ± 0.0099 


8.38 ±5.30 


1.47 ±1.09 


1.48 ± 1.10 


AEIo 


9/11 


0.0155 ± 0.0098 


11.38 ±6.18 


6.55 ±2.73 




FAU 


9/11 


0.0021 ± 0.0035 


8.51 ±4.75 


3.67 ± 1.68 




BKL 


9/11 


0.0153 ± 0.0094 


11.74 ±6.39 


6.89 ±2.92 




Kes 


9/11 


0.0174 ±0.0105 


11.99 ±6.51 


7.04 ± 2.96 




AEIn 


2/3 


0.0205 ±0.0127 


11.96 ±6.17 


1.81 ± 1.21 


1.83 ± 1.24 


AEIo 


2/3 


0.0199 ±0.0124 


14.10 ±6.99 


7.37 ±2.80 




FAU 


2/3 


0.0034 ± 0.0026 


10.66 ±5.61 


4.41 ± 1.76 




BKL 


2/3 


0.0191 ±0.0108 


14.52 ± 7.05 


7.79 ± 2.98 




Kes 


2/3 


0.0217 ±0.0124 


14.83 ± 7.24 


8.02 ±2.99 




AEIn 


1/3 


0.0165 ±0.0109 


8.58 ±4.17 


3.96 ± 4.46 


4.25 ±5.16 


AEIo 


1/3 


0.0156 ±0.0101 


9.57 ±4.62 


10.45 ±4.12 




FAU 


1/3 


0.0177 ±0.0090 


7.80 ±3.84 


6.60 ±2.75 




BKL 


1/3 


0.0148 ± 0.0089 


9.81 ±4.79 


11.24 ± 4.45 




Kes 


1/3 


0.0167 ±0.0105 


10.01 ±5.06 


11.49 ±4.48 





